Draft version February 4, 2008 

Preprint typeset using I4TgX style emulateapj v. 6/22/04 



MICROLENSING SENSITIVITY TO EARTH-MASS PLANETS IN THE HABITABLE ZONE 

Byeong-Gon Park', Young-Beom Jeon', Chung-UkLee' and Cheongho Han^-^ 

Draft version February 4, 2008 

ABSTRACT 

Microlensing is one of the most powerful methods that can detect extrasolar planets and a future space- 
based survey with a high monitoring frequency is proposed to detect a large sample of Earth-mass planets. In 
this paper, we examine the sensitivity of the future microlensing survey to Earth-mass planets located in the 
habitable zone. For this, we estimate the fraction of Earth-mass planets that will be located in the habitable zone 
of their parent stars by carrying out detailed simulation of microlensing events based on standard models of the 
physical and dynamic distributions and the mass function of Galactic matter From this investigation, we find 
that among the total detectable Earth-mass planets from the survey, those located in the habitable zone would 
comprise less than 1% even under a less-conservative definition of the habitable zone. We find the main reason 
for the low sensitivity is that the projected star-planet separation at which the microlensing planet detection 
efficiency becomes maximum (lensing zone) is in most cases substantially larger than the median value of the 
habitable zone. We find that the ratio of the median radius of the habitable zone to the mean radius of the 
lensing zone is roughly expressed as c/hz/'"e ^ 0.2(ot/0.5 Mq)'/^. 

Subject headings: planetary systems - planets and satellites: general - gravitational lensing 



1. INTRODUCTION 

Since the first detection by using the pulsar tim- 
ing method (^olszczan & Frail 1992), extrasolar plan- 
ets have been and are going to be detected by us- 
ing various techniques including the radial velocity tech- 
nique (Mayor & Oueloz 1995; Marcv & Butler 1996), transi t 
method (Struve 1952), astrometric technique (Sozzetti 2005'), 
direct imaging ( Anael 1994; Stahl & Sandler 1995), and mi- 
crolensing (Mao & Paczvhski 1991; Gould & Loebl 119921) . 
See the reviews of Perrvman (2000, 2005). 

The microlensing signal of a planet is a short-duration per- 
turbation to the smooth standard light curve of the primary- 
induced lensing event occurring on a background source star. 
The planetary lensing signal induced by a giant planet with 
a mass equivalent to that of the Jupiter lasts for a duration 
of 1 day, and the duration decreases in proportion to the 
square root of the mass of a planet, reaching several hours for 
an Earth-mass planet. Once the signal is detected and ana- 
lyzed, it is possible to determine the planet/primary mass ra- 
tio, q, and the projected primary-planet separation, s, in units 
of the Einstein ring radius te, which is related to the mass of 
the lens, m, and distances to the lens and source, Di and Ds, 
by 

Currently, several experiments are going underway to search 
for planets by using the microlensing technique and two ro- 
bust detections of Jupiter-mass planets were recently reported 
by Bond et al. (2004) and Udalski et al. (2005). In addition, 
a future space-based survey with the capability of continu- 
ously monitoring stars at high cadence by using very large 
format imaging camer as is propose d to detect a large sample 
of Earth-mass planets dBennett & Rhie.2002) . 
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The microlensing technique has various advantages over 
other methods. First, microlensing is sensitive to lower-mass 
planets than most other methods and it is possible, in prin- 
ciple, to detect Earth-mass planets from ground-based obser- 
vations wi th adequate monitoring fr equency and photometric 
precision ('Goul d. Gaudi. & Han ^2004). Second, a large sam- 
ple of planets, especially low-mass terrestrial planets, will 
be detected at high S/N with the implementation of future 
lensing surveys, and thus the microlensing technique will 
be able to provide the best statistics of Galactic population 
of planets. Third, the planetary lensing signal can be pro- 
duced by the planet itself, and thus the microlensing tech- 
nique is the only proposed method that can detect a nd charac- 
terize free-floating planets (bBennett & Rhie, .20021: iHan et all 
l2005h . Fourth, the microlensing technique is distinguished 
from other techniques in the sense that the planets to which 
it is sensitive are much more distant than those found with 
other techniques and the method can be extended to search 
for planets located even in other galaxies ( ICovone et alJl2000: 
Baltz & Gondolo 2001). 

In addition to the detections of planets, the habitability of 
the detected planets is of great interest. One basic condition 
for the habitability is that planets should be located at a suit- 
able distance from their host stars. Then, a question is whether 
a significant fraction of planets to be detected by future lens- 
ing surveys would be in the lensing zone. If the planets in 
the habitable zone comprise a significant fraction of the total 
microlensing sample of planets, a statistical analysis on the 
frequency of terrestrial planets in the habitable zone would 
be possible under some sorts of assumptions about the shape 
and inclination of the planetary orbit around host stars. In this 
paper, we examine the sensitivity of future lensing surveys 
to Earth-mass planets located in the habitable zone by car- 
rying out detailed simulation of Galactic microlensing events 
based on the standard models of the physical and dynamic 
distributions and the mass function of Galactic matter The 
sensitivities to habitable-zo ne planets of other planet search 
methods were discussed by Gould. Pepper. & DePov (2003|) 
(transit meth od), Gould, Ford, & F ischer (2003) (astrometric 
method), and lSozzetti et alJ (l2002h (space interferometry). 
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The paper is organized as follows. In § 2, we briefly de- 
scribe the basics of planetary microlensing. In § 3, we de- 
scribe the procedure of the simulation of Galactic microlens- 
ing events produced by lenses with Earth-mass planet com- 
panions. We then describe the procedure of estimating the 
fraction of planets located in the habitable zone among the 
total number of detectable Earth-mass planets based on the 
planetary lensing events produced by the simulation. In § 4, 
we present the results of the simulation and discuss about the 
results. We conclude in § 5. 

2. BASICS OF PLANETARY MICROLENSING 

Planetary lensing is an extreme case of binary lensing with 
a very low-mass companion. Because of the very small mass 
ratio, planetary lensing behavior is well described by that of 
a single lens of the primary for most of the event duration. 
However, a short-duration perturbation can occur when the 
source star passes the region around caustics. 

The caustic is an important feature of binary lensing and it 
represents the source position at which the lensing magnifica- 
tion of a point source becomes infinite. The caustics of binary 
lensing form a single or multiple closed figures where each of 
which is composed of concave curves (fold caustics) that meet 
at cusps. For a planetary case, there exist two sets of discon- 
nected caustics. One small 'central' caustic is located close to 
the primary lens, while the other bigger 'planetary' caustic(s) 
is (are) located away from the primary at the position with a 
separation vector from the primary lens of 

r = s(l-i) , (2) 

where s is the position vector of the planet from the primary 
lens. Then, the caustics are located within the Einstein ring 
when the planet is located within the separation range of 0.6 < 
s < 1.6, which is often referred as the 'lensing zone'. The 
number of the planetary caustics is one or two depending on 
whether the planet lies outside (s > 1) or inside (s < 1) the 
Einstein ring. 

The size of the caustic, which is directly proportional to 
the planet detection efficiency, is dependent on both q and 
s. Under perturbative approximation (when ^ <C 1 and \s- 
1 1 ^ ^), the sizes of the central (Arcc) and planetary (Arpc) 
caustics as measured by the width along the star-planet axis 
are represented, respectively, by 

Arcc 4g (4qs-^, for,s>l, 

re ~ (i-i-i)2 ^ \4qs\ fori < 1, ^ ' 

and 

Arpc ^ Uq^/^[s(s^- 1)'/^' -^4^'''V2, for,? > 1, 

rE ~\2g''/2(Ko-Ko' + Ko.r2)cos6lo^ L3^'/V, for5<l, 

(4) 

where ^(6*) = { [cos 26 ± (/ - sin^ 261)'/ ^] /(s^ - 1 /jf^)} '/^ Op = 
j2r±siirW3£/2)l /2, and kq = k(6Io) (' Bo77al2nn(lfAnl2nn.'Tl: 
lChungetal.1 12005: Han 2006). The size of the caustic be- 
comes maximum when i ~ 1 and decreases rapidly as the sep- 
aration becomes bigger (cx for both the central and plane- 
tary caustics) or smaller (cx for the central caustic and oc 
for the planetary caustic) than the Einstein ring radius. As a 
result, only planets located within the lensing zone have non- 
negligible chance of producing planetary signals. 

In Figure [2 we present the variation of the total caustic 
size. Arc = Arcc + Arpc, as a function of the primary -planet 




Fig. 1 . — The variation of the caustic size as a function of the primary- 
planet separation, s, for a planetary lens system with a mass ratio of g = 10"'. 
The dotted curve is calculated numerically while the dashed curve is com- 
puted by usingthe analytic formalism based on perturbative approximation 
(Eqs. (51 and 0). The two results match very well except the region around 
.V = 1, where the perturbative approximation is no longer valid. The solid 
curve is the model adopted in our simulation as the relative planet detection 
efficiency as a function of i. 



separation for a planetary lens system with a mass ratio of 
q = 10"^, which corresponds to an Earth-mass planet around a 
star with 0.3 Mq. In the figure, the dotted curve is calculated 
numerically while the dashed curve is computed analytically 
by using the formalism in equations Q and (0). The disagree- 
ment between the two curves in the region around = 1 is due 
to the failure of the perturbative approximation in this region. 

3. SEVIULATION 

To estimate the fraction of planets located in the habitable 
zone among the total number of Earth-mass planets detectable 
by future lensing surveys, we conduct detailed simulation of 
Galactic microlensing events. In the simulation, we assume 
that the survey is conducted toward the Baade's Window field 
centered at the Galactic coordinates of il,b) = (l°,-4°). The 
absolute brightnesses of the source stars are assigned on the 
basis of the luminosity function of |Holtzman et al. (199^ 
constructed by using the Hubble Space Telescope. Once the 
absolute magnitude is assigned, the apparent magnitude is de- 
termined considering the distance to the source star and ex- 
tinction. The extinction is determined such that the source 
star flux decreases exponentially with the increase of the dust 
column density. The dust column density is computed on the 
basis of an exponential dust distribution model with a scale 
height of h~ = 120 pc, i.e. cx exp{-\z\/hz), where z is the dis- 
tance from the Galactic plane. We normalize the amount of 
extinction so that A v = 1 .28 for a star located at Z).s = 8 kpc fol- 
lowing the measurement of lHoltzman etaTUIggl . The pro- 
posed space-based lensing survey searching for Earth-mass 
planets will monitor main-sequence stars to minimize the 
finite-source effect, which washes out the planetary lensing 
signal. We, therefore, assume that the brightness range of 
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source stars to be monitored by the survey is 18.5 < V < 25.0, 
which corresponds to early F to early M-type main-sequence 
stars. 

The locations of the source stars and lens matter are al- 
located based on the standard mass distribution model of 
|Han & Gould (2003). In the model, the bulge mass distribu- 
tion is scaled by the deprojected infrared light density profile 
of iDwek et aL (J995i) . specifically model G2 with R^^^ = 5 
kpc from their Table 2. The velocity distribution of the 
bulge is deduced from the tensor virial theorem and the re- 
sulting distribution of the lens-source transverse velocity is 
listed in Table 1 of Han & Gould (1995), specifically non- 
rotating barred bulge model. The disk matter distribution 
is modeled by a double-exponential law, which is expressed 
as p(R,z) = /5oexp[-(r-/?o)AR+ kl/Z^z], where (R,z) is the 
Galactocentric cylindrical coordinates, /?o = 8 kpc is the dis- 
tance of the sun from the Galactic center, po = 0.06 Mq pc~^ 
is the mass density in the solar neighborhood, and /is = 3.5 
kpc and /z. = 325 pc are the radial and vertical scale heights. 
According to the models of the dynamical and physical distri- 
butions, the ratio between the rates of bulge-bulge and disk- 
bulge events is Tb : Td = 61.8 : 38.2. 

We assign the lens mass based on the model mass func- 
tion of Gould (2000). The model mass function is composed 
of stars, brown dwarfs (BDs), and stellar remnants of white 
dwarfs (WDs), neutron stars (NSs), and black holes (BHs). 
The model is constructed under the assumption that bulge 
stars formed initially according to a double power-law dis- 
tribution of the form 



dm V m, 



-2.0 
-1.3 



for m > 
for m < nic 



(5) 



where nic = 0.7 Mr. These slopes are consistent with the 
observations of Zoccali et al. (2000) except that the profile 
is extended to a BD cutoff of m = 0.03 Mq. Based on 
this initial mass function, remnants are modeled by assum- 
ing that the stars with initial masses 1 Mq < m < 8 Mq, 
8 Mq < m < 40 Mq, and m > 40 Mq have evolved into 
WDs (with a mass 0.6 Mq), NSs (with a mass 13.5 Mq), 
and BHs (with a mass 5 Mq), respectively. Note that since 
we assume stars with m > 1 Mq have evolved into remnants, 
the upper limit of the stellar lens mass is 1 Mq, correspond- 
ing to the turn-off star in the Galactic bulge. ^ Then, the re- 
sulting mass fractions of the individual lens components are 
stars : BD : WD : NS : BH = 62 : 7 : 22 : 6 : 3. An important 
fraction of lenses are stars and thus they will contribute to the 
apparent brightness of the source star We consider this by de- 
termining the lens brightness by using the mass-My relation 
Usted in Allen (2000). 

For each event involved with a primary lens, we then in- 
troduce a companion of an Earth-mass planet. The loca- 
tion of the planet around the primary star is allocated under 
the assumption of a circular orbit with a power-law distri- 
bution of semi-major axes a, i.e. dNp/da ex. a~", and ran- 
dom orientation of the orbital plane. There is little con- 
sensus about the power of the semi-major axis distribution. 
IJabachnik & Tremaine (2002) claimed a = I from the analy- 
sis of observed extrasolar planets detected by radial velocity 

We note that disk turn-off is a bit brigliter. This may affect our result 
because it is easier to detect planets in the habitable zone with microlensing 
for more massive primaries. However, the effect will be small because (1) 
the minority of events are due to disk lenses, (2) the disk turn-off is not much 
brighter than the bulge one, and so (3) the number of potential main-sequence 
lenses with mass significantly higher than a solar mass is quite small. 



surveys. The minimum mass solar nebula has surface den- 
sity with a = 1 .5 (Hayashi 1995). Kuchner (2004) argued that 
multi -planet extrasolar planetary systems indicate a = 2.0. 
We, therefore, test three different powers of a= 1.0, 1.5, and 
2.0. The range of the semi-major axis is -2 < log(fl/AU) < 2, 
which is wide enough to cover the lensing zone of all possible 
lenses, and thus the exact range we consider does not affect 
the relative frequencies we calculate here. Once a planetary 
event is produced, the rate of each event is computed by 



Tpj p(Ds)Dlp(DOaLvAr„ 



(6) 



where p(D) is the matter density along the line of sight, the 
factor Dg is included to account for the increase of the number 
of source stars with the increase of Ds, ai represents the lens- 
ing cross-section corresponding to the diameter of the Ein- 
stein ring, i.e. ai = 2r^, and v is the transverse speed of the 
lens with respect to the source star The factor Arc is included 
to weight the planet detection efficiency by the cross-section 
of the planetary perturbation under the assumption that the 
planet detection efficiency is proportional to the caustic size.^ 
For the acceleration of the computation, we use Arc com- 
puted analytically by using the formalism derived under the 
perturbative approximation (Eqs. Js) and 0) except the re- 
gion where the approximation is not valid. In this region, we 
use the mean value of the numerical results averaged over the 
region (the solid curve in Fig.[0.^ Then, the fraction of events 
with planets in the habitable zone among the total number of 
events with detectable Earth-mass planets is calculated by 



/hz - 



I] rp.,(<iin < g < ^out) 

Erp,,- 



(7) 



where t/in and doui are the inner and outer limits of the habit- 
able zone, respectively. The semi-major axis is related to the 
projected separation by 



(cos^ (/)-!- sin 0cos^O'''^ 



(8) 



where is the phase angle of the planet measured from the 
major axis of the apparent orbit of the planet around the host 
star and / is the inclination angle of the orbital plane. 

The habitable zone is defined as a shell region around a star 
within which a planet may contain liquid water on its surface. 
A convenient way of estimating the inner and outer limits of 
the habitable zone is to assume that the planet behaves as a 
graybody with an albedo and with perfect heat conductivity 
(implying that the temperature is uniform over the planet's 
surface). Under this assumption, the inner and outer limits of 
the habitable zone are determined as 



.*in{out) ■ 



{1-Al)U 



1/2 



(9) 



' According to the formalism in equation the size of the caustic scales 
as .V"- for wide- separation planets {s 2> 1), whereas the detection probabil- 
ity actually scales as .v"'. Then, the assumption that the planet detection 
efficiency is proportional to the caustic size fails to apply to these wide- 
separation planets. However, we note that the habitable zone is, in most 
cases, much smaller than the Einstein ring radius, and thus this failure does 
not affect our result. 

^ The caustic size presented in Fig.QJis for an Earth-mass planet around 
a primary with a fixed mass of 0.3 Mq, while the events produced by the 
simulation are associated with primary stars of various masses. We note, 
however, that for a given mass of the planet, nip, the physical size of the 
caustic does not depend on the mass of the primary, nu, because Arc ~ 

Arpc oc q^f^rE oc (mflm^)^l^nij^ = nip. 
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Fig. 2. — The ranges of the habitable zone and lensing zone as a function 
of a stellar mass. For the habitable zone, the dark and light-shaded regions 
represent the conservative and less-conservative ranges of the habitable zone, 
respectively. The lensing zone is computed assuming that the distances to the 
lens and source star are Dl = 6 kpc and Ds = 8 kpc, respectively. 



where is the luminosity of the star, Al is the albedo of 
the planet, a is the Stefan-Boltzmann constant, and T is the 
radiative equilibrium temperature. In the simulation, we test 
two ranges of the habitable zone. For a conservative range, 
we adopt radiative equilibrium temperatures of T = 269 K 
and 203 K for the inner and outer Umits of the habitable 
zone, respectively. A less-conservative range has the same 
inner radius as that of the conservative range, but the outer 
limit is defined by a lower radiative equ ilibrium temperature 
of r = 169 K ( Lopez. Schneider & Danchi 200 5). For the 
albedo, we choose Al = 0.2 as a representative value of an 
Earth-like planet. For a sun-like star, then, the inner limit 
of the habitable zone is = 0.96 AU and the outer lim- 
its are t/out = 1-69 AU and 2.43 AU for the conservative and 
less-conservative ranges, respectively. For the conversion of 
the mass of the sta r into the lu minosity, we use the mass- 
luminosity relation jAllenl2000 ') of 

log(A)=3.81og(£).0.08. (10) 

Figure |2 shows the conservative (dark-shaded region) and 
less-conservative (light-shaded region) ranges of the habitable 
zone as a function of the stellar mass. 

The habitable zone is defined only for planets associated 
with stellar population of lenses. Therefore, identifying the 
lens as a star is important to further isolate the sample of can- 
didate planets in the habitable zone. Fortunately, the proposed 
space-based lensing survey can sort out events associated with 
stellar lenses (hereafter bright-lens events) by analyzing th e 
flux from the lens star (iBennett & Rhi3l2()02t lHmll20()5l) . 
We, therefore, estimate an additional fraction of planets in the 
habitable zone among the planets involved with bright-lens 
events, /hz,*- We assume that a lens can be identified as a star 
if it contributes > 10% of the total observed flux. Under this 



Fig. 3. — The distributions of the semi-major axes of Earth-mass planets to 
be detected in future lensing surveys. The intrinsic semi-major axes of planets 
are assumed to be distributed by a power-law of the form dNp /da on a^°' . The 
thick curves with black (for a = 1.0), blue {a = 1.5), and red (a = 1.0) colors 
represent the distributions of all planets, while the curves with grey (a = 1.0), 
light-blue (a = 1.5), and pink (q = 2.0) colors are the distributions of planets 
involved with bright lenses. The curves shaded with the corresponding colors 
are the distributions of planets located in the habitable zone. 



condition, we find that bright-lens events account for ^22% 
of the total events (~ 17% of bulge-bulge and ^ 29% of disk- 
bulge events). Identifying the stellar nature of the lens is also 
important to determine the distance to the lens and mass of 
the lens. Once Dl and m are known, the projected star-planet 
separation in physical units is determined hy d = srg- 

4. RESULTS AND DISCUSSION 

In Table [3 we summarize the resulting fractions of events 
with planets located in the habitable zone determined from 
the simulation. In Figure |3l we also present the semi-major 
axis distributions of the Earth-mass planets to be detected by 
the future lensing survey. In the figure, the thick curves with 
black (for the power a = 1 .0 of the intrinsic semi-major axis 
distribution), blue {a = 1.5), and red {a = 2.0) colors repre- 
sent the distributions of all planets, while the curves with grey 
(a = 1.0), light-blue {a = 1.5), and pink (a = 2.0) colors are 
the distributions for planets involved with bright lenses. The 
curves shaded with the corresponding colors are the distribu- 
tions of planets located in the habitable zone. Here we use the 
less-conservative range. 

From the table and figure, we find the following trends. 

1 . Planets located in the habitable zone comprises a very 
minor fraction of all planets. We find that the fraction 
ranges /hz ^ 0.2% -1.3% depending on the model dis- 
tributions of the intrinsic semi- maior axes and the def- 
initions of the habitable zone. IBennett & Rhi3 ( l200l 
predicted that the total number of Earth-mass planets 
that can be detected from the proposed space-based 
mission is ^ 200 if every lens star has an Earth-mass 
planet at an optimal position of 2.5 AU. Even under 
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TABLE 1 

Fraction of Planets in the Habitable Zone 



disti'ibution definition of the planet fraction in the habitable zone 

model habitable zone out of total events out of bright-lens events 



dNp/daoca-^'^ conservative 0.16% 0.71% 

less-conservative 0.62% 2.70% 

dNp/daoca^'-^ conservative 0.26% 1.15% 

less-conservative 0.88% 3.80% 

dNp/da oca-^ ^ conservative 0.42% 1.74% 

less-conservative 1.29% 5.40% 



Note. — The fraction of habitable-zone Earth-mass planets to be detected in future 
lensing surveys. We present two fractions. One is the fraction out of the total number 
of planets (/hz) and the other is the fraction among the planets involved with bright- 
lens events for which the primary lenses can be identified as stars (/hz.*). For the 
definitions of the conservative and less-conservative ranges of the habitable zone, see 
§ 3. The models in the first column indicate the assumed distributions of the intrinsic 
semi-major axis of planets. 



this optimistic expectation, then, the number of plan- 
ets in the habitable zone will be at most ^ 3, which is 
not adequate enough for the statistical analysis on the 
frequency of terrestrial planets in the habitable zone. 

2. Planets are biased toward smaller semi-major axis with 
the increase of the power a of the intrinsic semi-major 
axis distribution (c.f. the curves drawn with black [a = 
1.0], blue [a = 1.5], and red [a = 2.0] colors). How- 
ever, the dependence of the resulting semi-major axis 
distributions of lensing planets on a is weak. 

3. Planets of bright lenses are biased toward larger semi- 
major axis compared to the distribution of the total 
planets. This is because bright lenses tend to be heavier 
and thus have larger mean radius of the lensing zone, 
i.e. re. 

4. Most of habitable-zone planets are associated with 
bright lenses (c.f. the pairs of curves shaded with black- 
gray, blue-light blue, and red-pink colors). This is be- 
cause the gap between the lensing zone and habitable 
zone becomes narrower with the increase of the lens 
brightness. If the sample is confined only to planets in- 
volved with bright lens events, as a result, the fraction 
of habitable-zone planets increases into /hz * ^ 0.7% - 
5.4%. 

Then, why is the microlensing sensitivity to Earth-mass 
planets located in the habitable zone so low? First of all, 
the projected star-planet separation at which the planet de- 
tection sensitivity becomes maximum is in most cases sub- 
stantially larger than the median value of the habitable zone. 
This can be seen in Figure |2] where the ranges of the lensing 
zone and habitable zone are plotted as a function of the lens 
mass. Roughly, the median radius of the habitable zone is 
linearly proportional to the lens mass, while the mean value 
of the lensing zone is proportional to the square root of the 



mass. We find that the ratio of the median radius of the hab- 
itable zone to the mean radius of the lensing zone is roughly 
expressed as 



0.2 



1/2 



(11) 



rE VO.5M0, 
As a result, the two zones overlap only for stars more mas- 
sive than the sun (Di Stefano 1999). In addition, the habitable 
zone is defined by the intrinsic separation, while the lensing 
zone is defined by the projected separation. Considering the 
projection effect, then, the gap between the lensing zone and 
habitable zone further increases. 

5. CONCLUSION 

We investigated the microlensing sensitivity of future lens- 
ing surveys to Earth-mass planets located in the habitable 
zone by conducting detailed simulation of Galactic microlens- 
ing events. From the investigation, we found that the planets 
located in the habitable zone would comprise a very minor 
fraction and thus statistical analysis on the frequency of terres- 
trial planets in the habitable zone based on the sample would 
be very difficult. We find the main reason for the low sensi- 
tivity is that the projected star-planet separation at which the 
microlensing planet detection efficiency becomes maximum 
is in most cases substantially larger than the median value of 
the habitable zone. 
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